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Th=  arguments  Shat  leal  to  the  x^3  dependence  of  vote  *t 
x-2/3  dependence  of  axial  velocity  on  axial  distance  la  the  f 
ante  are  revleved.  Assuming  similarity  of  mean  velocity  prol 
dominant  terms  la  the  Havicr-Stctes  equations  for  turbulent  > 
used  to  derive  these  profiles  «ith  the  help  of  Bousslneso.-s  . 
hypotheaes.  The  resulting  profile  has  the  same  form  as  In  0 
taovn  tvo-dimensional  case.  It  Is  compared  vlth  an  assumed  * 
file  and  both  are  compared  tilth  the  Gaussian  profile  that  re 
corresponding  laminar  case.  The  form  of  the  eddy  viscosity 
tained  in  the  tvo  canes:  for  the  former  It  is  aero  on  the  a 
to  a  maximum  before  falling  to  aero  again  cn  the  ante  bound. 
for  a  Gaussian  profile  it  is  constant  across  a  cross  sectio. 
hulent  vale.  In  addition,  the  eddy  viscosity  has  on  X 
both  cases,  vhich  is  absent  In  the  tvo-dimensional  vote. 
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I.  Introduction 

The  trails  of  axl-symmetric  hypersonic  hlant  bodies  flying  through  the 
atmosphere  are  at  least  partially  turbulent.  Ve  expect  that  the  turbulence 
has  some  effects,  perhaps  pronounced,  on  their  radar  and  optical  images,  so  that 
ir  ue  aant  to  monitor  such  bodies,  It  becomes  important  for  hypersonic  flo»  theory 

to  he  able  to  predict  these  effects. 

At  the  present,  theory  has  very  little  to  say  about  these  turbulence  effects, 
ve  have  to  consider  the  non-homogeneous  shear  flow  of  a  compressible,  chemically 
reacting  gas.  Statistical  theories  of  turbulence  exist  only  for  homogeneous 
fluids,  and  even  then  there  are  feu  results  for  compressible  fluids.  On  the 
other  hand,  the  semi-emplrieul  and  heuristic  procedures  that  have  been  developed 
for  shear  flous  again,  or.  the  vhole,  avoid  the  compressible  ease,  nevertheless, 
since  these  seml-empirleal  procedures  are  all  ue  have.  It  la  the  purpose 
,resent  report  to  revleu  them  and  to  apply  the.  specifically  to  the  axi-symmctrlc 

vuXe.  This  is,  in  any  case,  a  necessary  first  step. 

Throughout  this  report  then,  ue  shall  assuma  that  ue  are  dealing  vith  n 

continuum  vlth  constant  density  and  constant  molecular  viscosity.  At  ever, 
stQ£e  ve  shall  compare  the  turbulent  and  the  laminar  cases.  Ve  begin  uith  an 
order  of  magnitude  analysis  of  the  terms  In  .be  Vavier-StoKes  equations.  Assum¬ 
ing  the  mean  flou  is  steady  and  Independent  of  the  azimuthal  angle  in  cylindrical 
coordinates  It  ulll  be  shoun  that  the  fluctuating  part  of  the  velocity. 
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U  *  a  -u*  ,  enters  the  equations  only  through  the  combination  vhere 

tL  bar  denotes  tine  averages  and  x  and  r  refer  to  the  axial  and  radial  direc¬ 
tions  respectively.  Next  ve  present  the  well-known  momentum  theorem,  which- 
relates  the  total  drag  on  the  body  to  an  integral  of  the  axial  velocity  over 
a  cross  section  of  the  wake,  and  of  which  extensive  use  is  made  subsequently. 

Some  elementary  considerations  then  lead  to  the  conclusion  thn*  t.ie  width  o 
the  turbulent  wake  increases  as  xl/3  and  the  velocity  on  the  axis  decreases  as 
x-2/3  (ln  the  laminar  case  these  dependencies  are  and  x  1).  Finally, 

assuming  similarity  of  velocity  profiles  in  the  far  wake  and  raking  use  of  the 
concepts  of  an  eddy  viscosity  and  of  Pranitl‘3  mixing  length,  ve  derive  the 

approximately  Gaussian  shape  of  the  mean  velocity  profiles. 

Recently  a  much  more  sophisticated  analysis  of  the  real  hypersonic  wake, 
taking  con./  resslbility  into  account,  has  been  made  by  Lees  and  Hromas^]. 

Lees  and  Hromas  make  extensive  use  of  the  concept  of  an  eddy  viscosity,  of 
similarity-type  agruments,  and  of  various  forms  of  the  momentum  theorem.  Their 
principal  innovation  is  to  treat  the  more  realistic  picture  of  an  inner  turbulent 
core  surrounded  by  an  outer  laminar  wake,  -which  in  turn  is  surrounded  by  the 
shock  shape  and  the  undisturbed  medium.  As  the  core  grows  it  eventually  fills 
all  of  the  wake,  so  that  asymptotically  their  picture  tends  to  the  picture  con¬ 
sidered  here.  Their  results  agree  well  with  the  shadowgraph  data  of  Slattery 

and  Clayj^J  . 

n.  order  Of  Magnitude  Analysis  of  the  Terms  of  the  Havier-jtokes  Equation 

The  Navler-Stokes  equations  and  the  continuity  equation  for  steady  flow 
with  constant  density  and  constant  molecular  viscosity  in  cylindrical  coordinates 

(x,  r,  ♦)  are, 
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For  laminar  a*l-sy=etrlc  flo»  ve  put  u,  =  0  and  W  -  0.  It  ve  also  replace 
ux  by  u  +  dx  vhere  U  la  the  constant  Incident  velocity  of  the  fluid  In  the  rest 
frame  of  the  body,  then  for  such  laminar  flov  (1)  to  (4)  become, 
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vhere  an  important  simplification  has  been  effected  by.  dropping  all  terms  involv¬ 
ing  the  molecular  viscosity  in  (l),  (2),  and  <3>-  This  is  the  first  of  many 
assumptions  arising  from  empirical  knowledge.  The  tenns  involving  mean  values  . 
of  the  fluctuating  velocity  components  in  (8),  (9)*  and  the  tuibulent 

stresses,”  are  known  to  be  several  orders  of  magnitude  greater  than  the  molec¬ 
ular  viscous  stresses.  The  data  of  Townsend  [h]  ,  for  instance,  on  the  flow 
behind  an  infinite  cylinder  show  that  the  mean  properties  of  the  wake  are  inde¬ 
pendent  of  Reynolds  number  for  Reynolds  numbers  R  =  ^  greater  than  about  800, 
d  being  the  diameter  of  the  cylinder.  This  justifies  the  neglect  of  the  molec¬ 
ular  viscous  stresses.  On  the  other  hand,  we  also  know  that  molecular  viscosity 
is  the  only  available  mechanism  for  the  final  decay  of  the  turbulence.  Any 
possible  exact  treatment  of  equations  (8)  to  (ll)  would  therefore  be  incapable 
of  describing  the  process  of  dissipation.  In  fact,  of  course,  equations  (8) 
to  (ll)  cannot  be  solved  without  making  some  semiempirical  assumptions  about 
the  form  of  the  turbulent  stresses,  and  these  assumptions  will  implicitly  con¬ 
tain  information  on  the  effects  of  the  molecular  viscous  stresses  as  well. 

The  analogy  between  the  molecular  viscous  stresses  in  equations  (l)  to  (3) 
and  the  turbulent  stresses  in  equations  (8)  to  (10)  is  most  clearly  seen  by 
writing  the  Navier-Stokes  equations  compactly  in  rectangular  coordinates.  Using 
the  summation  convention  and  letting  i  -  1,  2,  3  refer  to  x,  y,  z  coordinates, 
the  laminar  equations  are 


dul  du^ 
3T  *uj  XT. 


VlJ 


* 


.  ..v  •**■'<*. 


where 


3u.  5a. 

(t<Aj  m  “  (5^  +  3jJ> 


and  the  turbulent  equations  without  terms  involving  the  molecular  viscosity  are 


where 
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The  laminar  equations  form  a  closed  system:  three  equations  for  three  unknowns. 
The  turbulent  equations  involve  in  addition  to  the  six  independent  components 
of  the  turbulent  stress  tensor  (r^.)^.  If  we  hope  to  solve  the  latter  we  must 
make  some  assumptions  relating  (T^.)^  to  the  mean  velocity  or  its  gradient. 

To  effect  an  order  of  magnitude  analysis  of  equations  CO)  to  (ll)  we  fol¬ 
low  Hinze  ^2j  and  introduce  longitudinal  and  transverse  length  scales,  and 
I <r>  which  are  distances  characteristic  of  significant  changes  in  the  mean  flow. 
For  wake  flow,  expecially  at  high  velocities,  we  assume 
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If  Ux  and  are  the  corresponding  longitudinal  and  transverse  velocity  scales 
(Ux  gives  the  scale  for  u^  only,  i.e.  it  does  not  include  the  velocity  at  infin¬ 
ity  U),  then  it  follows  from  the  continuity  equation  (ll)  that 
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ut  /  the  scale  for  uf,  of,  Zf,  which  are  all  assumed  to  he  of  the  same 
order  of  magnitude  (c.f.  Tovneend  [k]  for  the  slightly  different  ease  of  an 


infinite 


cylinder).  If  correlation  coefficients  arc  defined  by 
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then  the  scale  site  of  uj  o’  .  In  general  .111  he  8^  V  . 

With  these  seel?  Bi7.es  the  orders  of  magnitude  of  the  terns  in  (3)  are 


(after  multiplying  through  by  LjV%  ), 
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and  those  in  (9)  are  (after  multiplying  through  by  Wx* 


(AP)  l  2  v2  Lx  R  rL  ,(15) 
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where 


(Ap)  and  (Ap)r  are  the  gross  longitudinal  and  lateral  pressure  differ- 


ences. 


We  no.  mate  the  far  wa*e  assumption,  which  Is  that  the  difference  in  axial 
velocity  on  the  axis  and  very  far  Iron  the  axis  is  small,  1-e. 
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The  data  of  Townsend  M  ,  again  for  an  infinite  cylinder,  show  that  this  leads 
to  no  contradiction  for  axial  distances  downstream  from  the  body  greater  than 
about  100  body  diameters. 

If  we  make  the  further  reasonable  assumptions  that  Rrx'—0(l),  and  that 
v2/Ux&  is  at  most  0(l),  then  it  follows  from  ( 14)  and  (12)  that  U/Ux  is  at  most 
of  order  L JLt.  Therefore,  in  (15),  L.H.S.«  R.H.S.,  and  therefore,  the  pres 
sure  term  on  „he  R.H.S.  of  (15)  must  balance  the  fir6t  turbulence  term,  (which 
represents  the  second  and  fourth  terms  on  the  R.H.S.  of  (9))*  Thus  (9)  becomes 
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which  upon  integration  gives 


P/p 


=  const. 


(17) 


Differentiating  now  with  respect  to  x  gives 
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which  shows  that  the  pressure  term  in  (i.4)  is  of  order  v  / Ux  ,  and  may  therefore 
be  neglected  in  comparison  with  the  second  turbulence  terms.  (l4)  then  shows 
that  (8)  becomes  (after  dividing  through  by  U2), 
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This  is  the  equation,  that  will  be  used.  in.  the  sequel  to  derive  the  mean  proper¬ 
ties  of  the  turbulent  wake. 

Eauation  (17)  would  give  information  on  the  mean  pressure  distribution  if 
suitable  assumptions  about  u^  could  be  introduced.  We  also  have  not  yet  con¬ 
sidered  equation  (10);  the  order  of  magnitude  analysis  of  (10)  shews  that 
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i.e.,  the  axial  component  of  the  turbulent  velocity  is  much  more  strongly  cor¬ 
related  with  the  azimuthal  component  than  is  the  radial  component.  This  seems 
to  suggest  that  eddies  with  their  axes  parallel  to  the  main  stream  direction 
are  less  prevalent  than  eddies  with  their  axes  pointing  in  a  radial  direction. 

If  we  make  a  similar  analysis  of  the  equations  for  the  laminar  wake,  (5), 
(6),  (7),  then  instead  of  (l8)  we  get 
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III.  The  Momentum  Theorem 

In  this  section  wc  closely  follow  the  treatment  of  Landau  and  Lifshitz  £l 
and  write  the  Navier-Stokcs  (laminar)  equations  (i),  (2),  (3)  in  rectangular 
coordinates  in  the  form 


5 

5t 


(p  UL) 


where 


is  the  momentum  flux  density  tensor. 
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.here-  the  integral  1.  over  a  cross  section  of  the  viscoun  wake.  But  inside 
the  vake  p~pl2<Sf  P%.  a"*  «»  vlacoua  term  la  also  small,  especially  at 
hieh  Reynolds  numbers  (ve  are  still  considering  a  laminar  vake,  hovever). 


Thus  the  drag  becomes 
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If  the  vake  is  turbulent'  then  all  velocities  In  (20)  must  he  replaced  by 
dean  velocities,  and  the  molecular  viscous  tern.  In  (20)  must  be  replaced  by 
v  ^  ,  n,e  analysis  then  gees  through  as  for  the  laminar  vake,  end  (22) 


becomes 
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VJ7)  shows  that  pU  ux  is  again  the  dominant  term,  eo  that 
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In  the  next  section  we  shall  make  immediate  use  of  (23)  and  (24), 


IV.  The  Width  of  the  Wake  and  the  Velocity  on  the  Axis 
A.  Laminar  Wrke 

Let  b(x)  be  the  width  of  the  vake,  i.e.  assume  that  at  any  fixed  x  u^ 
becomes  negligible  for  r  >b(x).  Then  in  the  equation  for  the  laminar 
vake,  (19),  the  order  of  magnitude  of  the  L.H.3.  Is  uj<.x-x„)  U,  where  xo 
is  some  virtual  origin  of  the  vake  (close  to  the  body).  The  order  of  mag 
ultude  of  the  B.H.S.  of  (19)  1.  ru^/bV.  Equating  the  tvo  gives, 

1/2  1/2 

b^(x)  (|j)  (x_x0) 


(25) 
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Further,  from  the  momentum  theorem  1^3)  we  hare 
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vhere  u  is  the  velocity  on  the  axis,  therefore 
xo 


u  *  I  1 

_*2.  —  — ^  (x-x  )  • 

j  pyU  o 


B.  Turbulent  Wake 

Before  any  information  about  the  turbulent  wake  can  be  extracted  from 
equation  (l8)  some  assumption  relating  the  turbulent  shear  stress  Purux 
to  the  mean  velocity  must  be  made.  The  assumption  that  is  usually  made 
is  implied  in  the  terminology  that  we  have  already  been  using:  in  analogy 
with  the  molecular  viscous  stresses  we  put, 
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This  is  Boussinesq’ s  hypothesis,  l  is  called  the  dynamic  eddy  (or  turbu¬ 
lent)  viscosity.  In  contrast  to  the  laminar  case  we  cannot,  in  general, 
go  on  to  assume  that  e  is  an  intrinsic  property  of  the  fluid,  even  for  an 
incompressible  fluid.  To  relate  the  eddy  viscosity  in  turn  to  the  mean 
properties  of  the  flow  the  usual  procedure  is  to  use  Prandtl’s  mixing 
length  hypothesis.  In  analogy  to  the  mean  free  path  of  gas  molecules  this 
introduces  the  concept  of  a  length  l,  characteristic  of  the  decay  of  fluid 
"particles”  (i.e.  eddies).* 

*  A  good  discussion  of  Prandtl’s  mixing  length  hypothesis  is  given  by 
Schllchting  (3)»  P*  ^77. 
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(28] 


and  we  make  the  further  reasonable  assumption  that  the  mixing  length,  l 
is  of  the  same  order  as  the  width,  t>t,  of  the  turbulent  wake. 


With  the  three  assumptions  (27),  (28),  (29)  the  order  of  magnitude 
-of  the  R.1I.S.  of  equation  (l8)  becomeB 


and  equating  this  to  the  order  of  magnitude  of  the  L.H.S.,  \iJ(x-xq)  Uj 
give  8 
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From  the  momentum  theorem  (24)  we  have 
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Therefore  ( 30 )  an<l  (3l)  finally  give 
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C.j.rw,  (32)  and  (33)  vith  (25)  and  (26)  notion  that  the  Oaalnar 
.aha  videns  .ore  rapidly  and  its  .ean  aria!  velocity  falls  off  .ore  rapidly 
than  the  turbulent  vahe.  me  may  see.  paradoxical,  since  one  vould  expect 
that  turbulent  dissipation,  Vhlch  is  so  much  greater  than  .olecular  viscous 
dissipation,  vould  produfce  a  shorter  turbulent  vahe.  Actually,  in  order 
to  define  the  "length"  of  a  vahe  one  vould  have  to  introduce  some  cut-off 
criterion,  such  as  .ensuring  the  vahe  only  to  the  paint  at  mch  the  vel¬ 
ocity  on  the  axis  has  reached  99#  of  the  incident  stream  velocity.  Hence 
the  relative  magnitudes  of  the  coefficient,  of  the  povers  of  X  shove  are 
important.  Unfortunately  the  theory  does  not  give  any  clues  as  to  these 

magnitudes. 

It  should  be  mentioned  that  Landam  and  Lifshitz  [l]  ,  P-  136,  derive 
(32)  and  (33)  by  different  reasoning  without  invoking  either  Boussinesq  0 
or  Prandtl's  hy^thesis.  Yet  another  approach  is  that  outlined  below. 

V.  MPim  Velocity  Distribution  in  the  Vake 

If  ve  assume  that  the  mean  flov  pattern  in  the  far  vahe  has  self-preserving 

form  ve  can  determine  vhat  that  form  is.  Here  ve  follov  the  procedure  of 

Hinte  [2]  ,  vho  treated  the  slightly  different  case  of  a  tvo  dimensional  vahe. 

Ms*,  a  transformation  of  caution  (IS)  from  the  independent  variable. 

(x,  r)  to  the  new  independent  variables  (|,  b)/  where 
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and  put 
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Thus  (l8)  "becomes 
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With  the  same  substitutions  the  momentum  theorem  (24)  gives 
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or 
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const . 
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(37) 


where 


(38) 


The  integral  in  the  momentum  theorem  is  over  an  infinite  plane  cross  section 

of  the  wake,  but  since  by  (35)  the  integrand  is  essentially  negligible  outside 

the  wake  r>b.  (x),  the  integration  can  be  taken  from  0  to  1  as  in  (38). 

x, 

Eliminate  t  Tran  (37)  and  (36)  so  that  (36)  may  be  written  in  the  fom 


2®t 
t  df 


(39) 


Since  the  L.H.S.  of  (39)  is  a  function  of  \  only,  and  the  R.H.S.  is  a  function 
of  11  only,  both  sides  must  equal  a  constant.  Let 
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